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SECTION  I 


INTRODUCTION  AND  COMPUTER  PROGRAM  DESCRIPTION 


This  report  is  the  sixth  and  last  volume  of  the  final  report  generated 
under  Air  Force  contract  F336 1 5— 80— C— 32 1 1  entitled  "Advanced  Life  Analysis 
Methods".  The  objective  of  this  contract  is  to  develop  the  design  criteria 
and  analytical  methods  necessary  to  ensure  the  damage  tolerance  of  aircraft 
attachment  lugs.  In  this  contract,  an  extensive  analytical  and  experi¬ 
mental  investigation  was  conducted  to  characterize  and  predict  fracture  and 
growth  behavior  of  cracks  in  attachment  lugs.  The  titles  of  the  six  volume 
final  report  generated  under  this  contract  are  listed  below. 


Volume  I  Cracking  Data  Survey  and  NDI  Assessment  for  Attachment 
Lugs. 

Volume  II  Crack  Growth  Analysis  Methods  for  Attachment  Lugs. 

Volume  III  Experimental  Evaluation  of  Crack  Growth  Analysis  Methods 

for  Attachment  Lugs. 

Volume  IV  Tabulated  Test  Data  for  Attachment  Lugs. 

Volume  V  Executive  Summary  ano  Damage  Tolerance  Criteria  Recommen¬ 
dations  for  Attachment  Lugs. 

Volume  VI  User's  Manual  for  "LUGRO"  Computer  Program  to  Predict 
Crack  Growth  in  Attachment  Lugs. 


This  volume  (Volume  VI)  describes  the  crack  growth  analysis  computer 
program  LUGRO  developed  under  this  contract  and  provides  user  input 
instructions  and  some  sample  input  and  output  data. 

LUGRO  is  an  automated  computer  program  which  has  been  developed  using 
the  state-of-the-art  methodologies  to  predict  the  residual  strength  and 
fatigue  crack  growth  behaviors  of  single  through-the-thickness  cracks  and 
single  corner  cracks  at  attachment  lugs.  ^his  section  further  describes 


the  operation  of  the  program,  functions  of  the  main  program  and  the 
subroutines,  and  a  flow-chart  of  the  computer  program. 

£i  1 ■'  '£>  yO.  /<£■ 
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The  program  includes  the  following  four  basic  elements: 
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o  Stress  intensity  factor  solution 

o  Baseline  crack  growth  rate  relationship 

o  Applied  load  sequence 

o  Spectrum  load-interaction  model 

Stress  intensity  factor  solutions  for  attachment  lugs  developed  under 

this  contract  for  a  variety  of  structural  and  loading  complexities  can  be 

found  in  Volume  II  of  this  final  report.  Some  of  these  solutions  have  been 

embedded  in  this  computer  program,  which  can  be  used  to  interpolate  the 

data  to  the  desired  lug  geometry.  The  program  also  has  the  option  to  input 

the  user-generated  stress/stress  intensity  data. 

The  optional  baseline  crack  growth  rate  relationships  embedded  in  the 

program  are  those  of  Paris,  Forman,  and  Walker.  The  optional  spectrum  load- 

interaction  models  incorporated  in  the  program  are:  the  Wheeler, 

Willenborg,  Generalized  Willenborg,  and  Hsu  models,  or  assuming  no  load 

interaction.  The  crack  growth  rate  equations  and  the  spectrum  load- 

interaction  models  are  discussed  in  detail  in  Volume  II  of  this  final 

report.  The  stress  levels  which  comprise  each  individual  mission  segment, 

and  from  which  a  mission  mix  spectrum  can  be  generated,  may  be  input  in 

five  optional  ways:  as  maximum  stress  (a  )  and  stresss  ratio  (R);  <t 

max  max 

and  minimum  stress  (cr  .  );  a  and  mean  stress  (<r  ) ;  a  and  alter- 

min  max  mean  mean 

nating  stress  (a  ^t);  or  R  and  stress  range  (Ao).  The  program  predicts 
the  crack  growth  using  a  block-by-block  integration  technique. 

For  tnrough-the-thickness  cracks,  either  the  compounding  solution  or 
the  Green's  function  solution  can  be  used  in  the  prediction.  In  predicting 
the  growth  behavior  of  a  single  corner  crack,  the  through-the-thickness 
crack  solution  may  be  modified  by  either  the  one-parameter  (compounding 
method  with  constant  a/c  ratio)  or  two-parameter  method  (Green's  function 
method).  For  one-parameter  analysis,  the  prediction  is  straightforward  and 
is  similar  to  through-the-thickness  crack  prediction.  For  two-parameter 
analysis,  it  is  assumed  that  for  a  given  number  of  applied  load  cycles,  the 
extension  of  the  quarter  elliptical  crack  border  is  controlled  by  the 
stress  intensity  factors  at  the  intersection  of  the  crack  periphery  at  the 
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hole  wall  and  the  lug  surface,  i.e.,  and  K^.  In  general,  the  stress 
intensity  factors  at  these  two  locations  are  different,  resulting  in 
different  crack  growth  rates.  Therefore,  the  new  flaw  shape  aspect  ratio 
after  each  crack  growth  increment  will  be  different  from  the  preceding  one. 
The  new  crack  aspect  ratio  is  computed  using  the  new  crack  lengths  on  both 
the  hole  wall  and  lug  surface.  The  process  is  repeated  until  the  crack 
length  along  the  hole  wall  is  equal  to  the  lug  thickness.  At  that  time  the 
transitional  crack  growth  criterion  is  used,  until  the  crack  has  achieved  a 
uniform  through-the-thickness  shape.  After  that,  if  the  failure  has  not 
occurred e  a  one-dimensional  through-the-thickness  crack  analysis  is  used  to 
continuously  predict  the  subsequent  crack  growth  life.  The  analysis  is 
considered  to  be  complete  when  fracture  occurs  due  to  fracture  toughness  or 
net-section  yield  criterion  or  when  the  desired  final  crack  length  or  the 
maximum  usage  time  is  reached. 

The  LUGRO  program  contains  one  main  program  and  14  subroutines.  The 
functions  of  the  main  program  and  the  subroutines  are  described  below. 

MAIN  -  The  main  program  receives  the  input  of  basic  data  for  the 

crack  growth  life  prediction;  interfaces  with  the  sub¬ 
routines;  and  outputs  the  predicted  crack  growth  history. 

DADNDK  -  This  subroutine  reads  the  baseline  constant  amplitude 
crack  growth  rate  data  in  the  form  of  da/dN  vs  Ak  or 
directly  as  material  constants.  In  the  former  case,  for 
each  straight  line  segment  of  da/dN  vs  4k  data  (in  log- 
log  scale),  the  subroutine  calculates  the  material  crack 
growth  rate  constants  using  the  crack  growth  rate  equa¬ 
tion  specified  by  the  user. 

SPECTM  -  This  subroutine  reads  and  generates  one  complete  block  of 
stress  spectrum  to  be  used  in  the  crack  growth  predic¬ 
tion.  This  block  of  stress  spectrum  will  be  repeated 
until  the  computed  maximum  stress  intensity  factor 
reaches  the  fracture  toughness  value  or  the  net-section 
stress  value  reaches  the  yield  strength  of  the  material. 


The  prediction  will  also  be  terminated  when  the  crack 
le"gth  or  the  accumulated  flight  time  in  terms  of  blocks 
reaches  the  specified  maximum  value  in  the  input. 

GROWA  -  This  subroutine  estimates  the  crack  length  as  a  function 
of  accumulated  flight  time  (or  usage).  The  computations 
for  through-the-thickness  cracks  are  made  in  this  sub¬ 
routine.  For  corner  cracks,  the  computations  are  made  in 
subroutine  CORNER,  which  is  called  by  this  subroutine. 
The  predictions  will  be  terminated  if  the  computed  K 

ilia  X 

reaches  Kc,  or  the  net-section  stress  value  reaches  the 
yield  strength  of  the  material,  or  the  current  crack 
length  or  the  accumulated  flight  time  exceeds  the  maximum 
value  specified  in  the  input. 

SIF  -  The  stress  intensity  factors  for  a  through-the-thickness 

crack  are  generated  in  this  subroutine  for  crack  lengths 
ranging  from  0.0001  inch  to  the  desired  maximum  crack 
length  specified  in  the  input.  This  subroutine  calls  the 
appropriate  subroutines  to  calculate  stress  intensity 
factors  using  the  compounding  method  or  Green's  function 
methods,  or  interpolating  the  tabulated  data  that  have 
been  input  or  are  already  available  internally  in  the 
program. 

SIF8  -  This  subroutine  computes  the  stress  intensity  factors  at 

the  intersection  points  of  a  quarter-elliptical  crack 
front  at  the  hole  wall  and  lug  surface  by  modifying  the 
through-the-thickness  solution  with  various  correction 
factors. 

SACKS  -  This  subroutine  estimates  a  back  surface  correction  fac¬ 
tor  for  a  part-through  corner  crack. 

CORNER  -  This  subroutine  calculates  the  two-dimensional  growth  of 
a  corner  crack  as  a  function  of  accumulated  flight  time 
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(or  usage).  This  subroutine  does  the  same  functions  as 
the  subroutine  GROWA,  except  that  this  is  for  a  corner 
crack.  Also  calls  the  appropriate  subroutines  and  com¬ 
putes  the  transition  and  subsequent  through-the-thickness 
crack  growths. 

COMPND  -  This  subroutine  computes  the  stress  intensity  factors 

using  the  simple  compounding  method. 

GREENF  -  This  subroutine  computes  the  stress  intensity  factors 

using  the  Green's  function  method. 

TRANSG  -  This  subroutine  computes  the  transitional  crack  growth 

from  when  the  corner  crack  breaks  through  the  back 

surface  to  when  the  crack  becomes  a  uniform  through-the- 
thickness  crack. 

DADP  -  This  subroutine  calculates  the  average  crack  growth  per 

block  of  stress  spectrum  as  a  function  of  any  given  crack 
length. 


DADNJJ  -  This  subroutine  calculates  the  crack  growth  rate  for  each 

computed  effective  stress  intensity  factor  range  ( <dKe^) 

and  ratio  (K  .  __/K  . 

mineff'  maxeff 

EFFKDA  -  This  subroutine  calculates  effective  stress  intensity 

factor  range  (  „.)  and  ratio  (K  .  -./K  and 

eff  mmeff  maxeff  ’ 

estimates  the  corresponding  crack  growth  increment  using 
the  specified  crack  growth  retardation  equation. 

DIRCAL  -  This  subroutine  contains  the  stress  intensity  factors  and 
the  stress  distributions  as  a  function  of  outer-to-inner 
radius  ratio  of  lug  in  a  tabular  form.  The  user  has  the 
option  to  use  these  tables  and  interpolate  the  data  for 
the  desired  lug  geometry. 


5 


Two  more  subroutines  are  referenced  in  the  program.  They  are  FUNCTION 
subroutines  GIRC  and  DTAB2,  which  are  UNIVAC  library  subroutines.  GIRC  is 
a  single  variable  Aitkin  interpolation  subroutine  and  DTAB2  is  a  two  vari¬ 
able  interpolation  subroutine. 

A  flow-chart  showing  the  interactions  between  the  main  program  and  the 
subroutines  is  shown  in  Figure  1-1.  The  main  program  first  processes  the 
material  crack  growth  data  in  DADNDK  and  develops  the  loading  spectrum  in 
SPECTM.  Then  it  starts  the  crack  growth  analysis  by  calling  the  subroutine 
GROWA.  The  through-the-thickness  stress  intensity  factor  solutions  are 
computed  first  by  calling  the  subroutine  SIF.  Based  on  the  input  specified 
by  the  user,  these  stress  intensity  factors  are  obtained  by  using  one  of 
the  various  methods  as  shown  in  the  flow-chart.  Then  depending  on  the 
crack  shape  (through-the-thickness  or  corner)  the  program  branches  to  com¬ 
pute  the  crack  growth  behavior  using  subroutines  DADP,  EFFKDA  and  DADNJJ. 
In  the  case  of  corner  cracks,  the  through-the-thickness  crack  solutions  are 
modified  with  various  correction  factors  for  corner  crack  stress  intensity 
factors  in  subroutine  SIF8.  Subroutine  TRANSG  computes  the  transitional 
crack  growth  behavior  of  corner  cracks. 
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0  USES  THE  UNIVAC  LIBRARY  FUNCTION  SUBROUTINE  GIRC 
@  USES  THE  UNIVAC  LIBRARY  FUNCTION  SUBROUTINE  DTAB2 


Figure  1-1.  Flow-Chart  of  the  Computer  Program  LUGRO 
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SECTION  II 


INPUT  INSTRUCTIONS 


This  section  describes  the  input  instructions  for  the  computer  program 
LUGRO.  The  program  can  be  used  to  predict  the  residual  strength  and 
fatigue  crack  growth  behaviors  of  single  through-the-thickness  cracks  and 
single  corner  cracks  at  attachment  lugs. 

The  program  is  capable  of  analyzing  attachment  lugs  with  and  without 
interference-fit  bushings.  Figure  2-1  shows  the  geometries  and  nomencla¬ 
tures  of  attachment  lugs  with  and  without  bushings.  In  the  case  of  corner 
cracks,  there  are  three  distinct  regions  of  crack  growth,  namely  corner 
crack  growth  until  the  crack  breaks  through  the  thickness,  transitional 
growth  until  the  crack  reaches  a  uniform  through-the-thickness  shape  and 
the  subsequent  through-the-thickness  crack  growth.  These  three  regions-  of 
crack  growth  and  corresponding  definitions  of  crack  length  are  illustrated 
in  Figure  2-2. 

The  input  instructions  for  the  computer  program  LUGRO  to  predict  the 
residual  strength  and  crack  growth  behaviors  are  provided  in  Table  2-1. 
There  are  25  different  lines  of  input  as  shown  in  the  table.  Some  lines 
may  have  to  be  repeated  several  times  depending  on  the  nature  of  input. 
These  are  indicated  under  the  column  ’COMMENT'.  For  example,  the  da/dN  vs 
4K  data  (Line  4)  have  to  be  repeated  'KPOINT'  times,  wnich  is  specified  in 
Line  3.  The  table  also  describes  the  variables,  type  of  variables  and  the 
columns  in  which  tne  data  are  input  for  all  the  25  lines  of  input.  Suffi¬ 
cient  instructions  are  provided  in  the  table  to  choose  the  right  set  of 
input  lines.  For  example,  appropriate  lines  3,  4  and  5  should  be  input 
based  on  the  value  of  'IEQU',  which  specifies  the  selected  crack  growth 
rate  equation.  The  table  also  describes  the  cases  when  there  will  be  no 
inputs  for  some  specific  lines.  For  example,  when  the  value  of  ’IBETA'  is 
equal  to  -1,  there  are  no  inputs  for  lines  20,  21  and  22. 

The  input  variables  specified  in  the  25  lines  of  Table  2-1  are 
described  in  detail  in  Table  2-2.  The  table  also  describes  some  of  the 
terminologies  used  in  the  description  and  two  important  notes  (Notes  A  and 
B)  at  the  end  of  the  table.  The  units  with  which  the  inputs  need  to  be 
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specified  are  also  given  in  the  table.  This  section  is  followed  by  the 
sample  problem  section  in  which  details  of  input  and  output  are  provided  to 
aid  the  user. 
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TRANSITIONAL  CRACK 


THROUGH-THE-THICKNESS 

CRACK 


Figure  2-2.  Definitions  of  Crack  Lengths 
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TABLE  2-1.  INPUT  INSTRUCTIONS  FOR  CRACK  GROWTH  ANALYSIS 
PROGRAM  'LUGRO'  FOR  ATTACHMENT  LUGS 


LINE 

COLUMN 

VARIABLE 

TYPE  OF  VARIABLE 

COMMENT 

i 

i  -SO 

TITLE 

ALPHANUMERICAL 

*** 

1-5 

6-10 

11-15 

16-20 

I ECU 

LDADN 

l.LC'AD 

LSEQ 

INTEGER- 

INTEGER 

INTEGER- 

INTEGER 

0,4,5 

<  A ) . 

IF 

IECiU  =  1, 

WALKER'S  EQUATION- 

IS  USED,  THEN 

1-5 

KPOINT 

INTEGER 

4 

1-10 

11-20 

EFFK ( I ) 
DADN(I) 

REAL 

REAL 

REPEAT  'KPOINT'  TIMES 

1-10 

1 1  -20 

CM 

RC 

REAL 

REAL 

w  >  4  ■)  z*  •> 

\B). 

IF 

I  ECU  =  2, 

FORMAN -'S  EQUATION 

IS  USED,  THEN 

**T; 

i  “J 
6-10 

KPOINT 

INCFNF 

INTEGER 

INTEGER 

4 

1-10 

11-20 

EFFK(I) 

DADN(I) 

REAL 

REAL 

INPUT  ONLY  IF  INCFNF  =  0, 
REPEAT  'KPOINT'  TIMES 

1-10 

11-20 

21-20 

FC(  I ) 

FN<  I  > 

EFFK ( I ) 

INPUT  ONLY  IF  INCFNF  =  i, 
REPEAT  'KPOINT-1 '  TIMES 

1-10 

11-20 

21-30 

FKC 

FR 

RC 

REAL 

REAL 

REAL 

3  *  4  ?  5 

*  c  > . 

IF 

I ECU  =  3, 

PARIS'  EQUATION  ) 

[S  USED,  THEN 

5 

i 

6-10 

KPOINT 

liMCPNP 

INTEGER 

INTEGER 

4 

1-10 
i  i  -20 

EFFK< I ) 
BABN ( I ) 

REAL 

REAL 

INPUT  ONLY  IF  INOPNP  =  0, 
REPEAT  'KPOINT'  TIMES 

4** 

1-10 

1 1-20 

2 1  -  30 

PC  (I) 

PN<  I) 

EFFK ( I ) 

REAL 

REAL 

REAL 

INPUT  ONLY  IF  INCPNP  =  1, 
REPEAT  '  KPOINT-1 '  TIMES 

c* 

1-10 

RC 

REAL 
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TABLE  2-1 


(CONTINUED) 


i-io 

1  1-20 
21-20 

THREKO 

CPITkC 

FTY 

REAL 

REAL 

REAL 

- 

1  -5 

I RETAR 

INTEGER 

6- 1 0 

I MODEL 

INTEGER 

0  INPUT  ONLY 

IF  IRETAR  =  1 

1  “*? 

I PLANE 

INTEGER 

.1—15 

EX  PM 

REAL 

EXPM  IS  INPUT  ONLY  IF  I MODEL  =  4 

Xf'MXTH 

REAL 

YkHXTH  AND  SOS  ARE  INPUT  ONLY 

£  OS 

REAL 

IF  I MODEL  =  2 

1-5 

M'3 

INTEGER 

4-10 

I  READ 

INTEGER 

11-13 

I LAYER 

INTEGER 

10 

1  -SO 

TITLEM 

ALPHANUMERICAL 

-  i 

1-10 

•3RATI0 

REAL 

11-20 

HOURM 

REAL 

^  * 

1-10 

SMAX 

REAL 

11-20 

R 

REAL 

INPUT  ONLY  IF  I READ® 1.  REPEAT 

21-30 

CPF 

REAL 

THIS  LINE  •'I LAYER'  TIMES 

12* 

1-10 

SMAX 

REAL 

11-20 

SMIN 

REAL 

INPUT  ONLY  IF  IR£AD=2,  REPEAT 

21-30 

CPF 

REAL 

THIS  LINE  /  I LAYER •"  TIMES 

11**^ 

1-10 

SMAX 

REAL 

11-20 

SMEAN 

REAL 

INPUT  ONLY  IF  I READ-3.  REPEAT 

21-30 

CPF 

REAL 

tHIS  LINE  ■'ILAYER-'  TIMES 

lO'1** 

1-10 

SMEAN 

REAL 

11-20 

SALT 

REAL 

INPUT  ONLY  IF  IREAD=4,  REPEAT 

21-30 

CPF 

REAL 

THIS  LINE  'ILAYER''  TIMES 

i 2**** 

1-10 

■3DEL 

REAL 

1 1  -20 

R 

REAL 

INPUT  ONLY  IF  IREAD=5»  REPEAT 

2 1  -20 

CPF 

REAL 

THIS  LINE  •'ILAYER-'  TIMES 

i  3 

REPEAT  LINES 

*'  TO  12  FOR 

ALL  POSSIBLE  MISSIONS,  AND  THEN 

INPUT  •'O'  ON 

COLUMN  5  TO 

SIGNAL  END  OF  MISSION  PROFILES. 

«  ,i 
i  -t 

1-5 

NPASSO 

INTEGER 

o—  i  0 

MSEOT 

INTEGER 

J  e- 

l  .» 

1-5 

NFLTS 

INTEGER 

4-10 

HIS I ON 

INTEGER 

REPEAT  THIS  LINE  •  NSEGT-'  TIMES 

1  4 
i  ; 

i  -5 

I  THRU 

INTEGER 

6-10 

METHOD 

INTEGER 

i  i-i? 

IBETA 

INTEGER 

Uz-20 

IOFTYP 

INTEGER 

I GFTYP  INPUT  ONLY  IF  METHOD  =  2 

OP  -2 
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TABLE  2-1.  (CONTINUED) 


,7 

1-10 

PhD  I  OS 

PEAL 

11-20 

WIDTH 

REAL 

21-20 

THICK 

REAL 

1 

1  -  i  o 

C RAC  KO 

REAL 

11-20 

:  RAC  f-.  F 

REAL 

1  •*> 

1-10 

A  020 I N 

REAL 

INPUT  ONLY  IF  I THRU  =  0 

22  INF 

UT  ONL /  IF 

I BETA  »  1 

1  -so 

TITLEB 

ALPHANUMERICAL 

- « 

1-5 

NBFT 

INTEGER 

B 

REAL 

1'? 

1-10 

COVERS 

REAL 

1 1-20 

SET  ATS 

REAL 

REPEAT  -'MBPT '  TIMES 

co  -  2i- 

TC 

I SETA  =  -1 

THERE  ARE  NO  INPUT 

3  FOR  LINES  20,  21-  22 

‘I* 5 » 

25  INPUT  ONLY  IF 

I BETA  =  0 

22 

1-5 

NO  10 

I NTEGER 

6-10 

I RS 1 0 

INTEGER 

11-15 

IWRITEA 

INTEGER 

24 

I- 10 

TSR 

REAL 

11-20 

EL 

REAL 

21-20 

EBEL 

REAL 

INPUT  ONLY  IF  IRSIG  =  -1  OR 

31-40 

DELS 

REAL 

METHOD  =  -2 

41-50 

POL 

REAL 

31-60 

FOB 

REAL 

25 

1-10 

YO 

REAL 

INPUT  ONLY  IF  IRSIG  =  0 

11-20 

XSIO 

REAL 

AND  METHOD  =  2 

REPEAT  'N5IG/  TIMES 

1-10 

YO 

REAL 

INPUT  ONLY  IF  IRSIG  =  1 

11-20 

XSIO 

REAL 

AND  METHOD  =  2 

2l  30 

RS 1 0 

PEAL 

REPEAT  NSIG'  TIMES 

25** 

1 1  -20 

.<510 

REAL 

INPUT  ONLY  IF  IRSIG  =  -1 

AND  METHOD  =  2 

REPEAT  -'ll  (-NSIG)  '  TIMES 

21-30 

RE  10 

REAL 

IMFUT  ONLY  IF  IRSIG  =  1 

AND  METHOD  =  -2 

REPEAT  '1  i  =  <NSIG)  TIMES; 

FOR  THIS  CASE,  IN  CARD  24 

ONLY  TBF:  AND  EBEL  NEED  TO  BE 
INPUT  USING  THE  SAME  FORMAT 

25*#** 

NO  LINE 

25  INRUT  IF 

IRSIG  =  -1  AND  ME 

THOD  =  -2 
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TABLE  2-2.  INPUT  VARIABLES  DESCRIPTION 


BEFORE  DESCRIBING  THE  INPUT  VARIABLES,  SOME  TERMINOLOGIES,  WHICH  ARE 
NOT  INPUT  VARIABLES,  ARE  DEFINED  HERE,  FIRST  TO  AID  THE  PREPARATION 
Or  LOAD  SPECTRUM  AND  TO-  CHOOSE  THE  OUTPUT  OPTION. 


BLOC!'  -  ONE  COMPLETE  PASS  T '■'ROUGH  THE  ENTIRE  SPECTRUM  BEFORE  ITS 

REPETITION.  ONE  BLOCK  CONSISTS  OF  MULT I -NUMBER  OF  SEGMENTS. 


SEGMENT  -  FRACTION  OF  A  BLOCK.  ONE  SEGMENT  CONSISTS  OF  ONE  FLIGHT  OR 
MULTIPLE  NUMBER  OF  CONSECUTIVE  FLIGHTS  OF  SAME  MISSION. 


MISSION  -  A  TYPE  OF  FLIGHT  IN  AIRPLANE  USAGE. 


FLIGHT  -  ONE  COMPLETE  PASS  THROUGH  A  GIVEN  MISSION.  ONE  FLIGHT 
CONSISTS  OF  MULT I -NUMBER  OF  LOAD  LAYERS. 


LAYER  -  ONE  ENTRY  IN  A  MISSION.  A  LAYER  CONSISTS  OF  A  NUMBER  OF 
CONSICUTIVE  CYCLES  OF  SAME  STRESS  AMPLITUDES. 


THE  FOLLOWING  DESCRIPTION  IDENTIFIES  ALL  OF  THE  INPUT  VARIABLES  AND 
THEIR  DEFINITIONS. 


TITLE  :  PROBLEM  IDENTIFICATION  TITLE. 


I ECU  :  OPTION  TO  USE  DIFFERENT  GROWTH  RATE  EQUATION. 
IF  TEQU  =  1,  WALKER'S  EQUATION  IS  USED, 

=  2,  FORMAN'S  EQUATION  IS  USED, 

=  3,  PARIS'  EQUATION  IS  USED. 


LDADN  :  OPTION  TO  PRINT  OUT  THE  INPUT  GROWTH  RATE  DATA. 

IF  LDADN  =  0,  PRINT  OUT  THE  GROWTH  RATE  DATA, 

=  1,  DO  NOT  PRINT  OUT  THE  GROWTH  RATE  DATA. 


l.I.OAD  :  OPTION  TO  ''PINT  OUT  THE  INPUT  LOADING  SPEC~RUM. 

IF  LLOAD  =  0,  PRINT  OUT  THE  INPUT  LOADING  SPECTRUM, 

=  1,  00  NOT  PRINT  OUT  THE  INPUT  LOADING  SPECTRUM. 


LrEQ  OPTION  TO  PRINT  OUT  THE  INPUT  FLIGHT  SEQUENCE  OF  A 

BLOCK  LOAD  SPECTRUM. 

IF  L E EG  =  «•'.  PRINT  00 ~  THE  INPUT  FLIGHT  SPECTRUM, 

=  I,  DO  NOT  Fc I NT  CUT  THE  INPUT  FLIGHT  SPECTRUM. 
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TABLE  2-2.  (CONTINUED) 


POINT 

EFFK  < I > 

DAE'N  <  i  i 

CM 

°C 

INCFNF 

FC  i  I ) 

FN(  I) 

FKC 

FR 

IMCFNP 

PC  (  I  ) 

PM  <  I  > 


NUMBER  OF  POINTS  ON  (DELTA  K  VS.  DA/DM)  CURVE •  BOTH  END 
POINTS  ARE  INCLUDED.  NOUMBER  OF  LINEAR  SEGMENTS  USED 
TO  REPRESENT  THE  GROWTH  RATE  CURVE  WILL  BE  =  KPOINT-1. 
NOTE  :  MAXIMUM  KPOINT  VALUE  IS  SO. 


DELTA  r.'  OR  EFFECTIVE  DELTA  K  CORRESPONDING  TO  EACH  INPUT 
POINT  <  UNIT  :  KSI*SQRT( INCH)  ) 


GROWTH  RATE  DA/DN  CORRESPONDING  TO  EACH  DELTA  K  OR  EFFECTIVE 
DELTA  K  (  UNIT  :  MICRO INCHES  ) 


CONSTANT  M'  USED  IN  WALKER'S  EQUATION. 


CUT-OFF  OF  the  MINIMUM  APPLIED  STRESS  RATIO.  I.E..  IF  RCRC,  R~P 


OPTION  TO  INPUT  EITHER  (  DELTA  K  b  DA/DN  )  PAIRS  OR  (  CF  b  NF  ) 
PAIRS  FOR  FORMAN'S  EQUATION. 

IF  INCFNF  =  0.  INPUT  (  DELTA  K  l>  DA/DN  )  PAIRS, 

=  1,  INPUT  (  CF  $<  NF  >  PAIRS. 

NOTE  :  Cr  ’  ?.(  'NF'  ARE  CONSTANTS  OF  FORMAN'S  EQUATION. 


CONSTANT  'C'  USED  IN  FORMAN'S  EQUATION. 
i  FOR  DELTA  K  IN  UNIT  KSI*SGRT ( INCH)  ) 


CONSTANT  N'  USED  IN  FORMAN'S  EQUATION. 
(  FOR  DELTA  K  IN  UNIT  KSI*SGRT< INCH)  ) 


KC  VALUE  TO  BE  USED  IN  THE  FORMAN'S  EQUATION. 
(  UNIT  :  KS I^SGRT ( INCH)  ) 


R-RATIO  OF  CONSTANT  AMPLITUDE  OF  INPUT  DA/DN  DATA. 


OPTION  TO  INPUT  EITHER  (  DELTA  K  b.  DA/DN  )  PAIRS  OR  (  CF  b  NF  ) 
FAIRS  FOR  PARTS'  EQUATION. 

IF  INCPNP  =  0,  INPUT  <  DELTA  K  b  DA/DN  )  PAIRS, 

=  1,  INPUT  <  CF  b  NF  )  PAIRS. 

NOTE  :  'CF'  b  'NF'  ARE  CONSTANTS  OF  PARIS'  EQUATION. 


CONSTANT  'C '  USED  IN  PARIS'  EQUATION. 

(  FOR  DELTA  K  IN  UNIT  KSI-*SQRT ( INCH )  ) 


CONSTANT  N-  USED  IN  PARTS'  EQUATION. 

'  FOR  DELTA  i  IN  UNIT  KSI*SGRT< INCH)  ) 
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TABLE  2-2.  (CONTINUED) 


THREKO  :  THRESHOLD  DELTA  K  OF  THE  MATERIAL  USED.  WHEN  CALCULATED 
DELTA  K  IS  LESS  THAN  OR  EQUAL  TO  THREKO,  NO  CRACK  GROWTH 
WILL  BE  CONSIDERED.  <  UNIT  :  KSI*SGRT < INCH)  ) 


CRITKC  s  CRITICAL  STRESS  INTENSITY  FACTOR  OR  FRACTURE  TOUGHNESS 
OF  THE  MATERIAL.  (  UNIT  :  KSI*SQRT( INCH)  ) 


FTY  :  TENSILE  YIELD  STRENGTH  OF  THE  MATERIAL.  <  UNIT  :  KSI  ) 


I RETAR  :  OPTION  TO  CONSIDER  SPECTRUM  RETARDATION  EFFECTS. 

IF  I RETAR  a  0,  NO  RETARDATION  EFFECT  WILL  BE  CONSIDERED, 

a  1,  SPECTRUM  RETARDATION  EFFECTS  WILL  BE  CONSIDERED. 


IMODEL  ;  OPTIONS  TO  USE  VARIOUS  RETARDATION  MODELS. 

IF  IMODEL  «  1,  WILLENBORG  MODEL  WILL  BE  USED, 

IF  IMODEL  =  2,  GENERALIZED  WILLENBORG  MODEL  WILL  BE  USED, 
*  3,  HSU  MODEL  WILL  BE  USED, 

=»  4,  WHEELER  MODEL  WILL  BE  USED. 


I PLANE  :  OPTIONS  TO  USE  PLANE  STRESS  OR  PLANE  STRAIN  IN  THE  CALCULATION 
OF  YIELD  ZONE. 

IF  I PLANE  *  1,  PLANE  STRESS  IS  USED, 

*  2,  PLANE  STRAIN  IS  USED. 


EX PM  :  WHEELER'S  RETARDATION:  MODEL  EXPONENT  M. 


XKMXTH  :  MAXIMUM  THRESHOLD  STRESS  INTENSITY  FACTOR  <  UNIT  *  KSI*SQRT( INCH)  ) 
(ONLY  FOR  GENERALIZED  WILLENBORG  MODEL) 


SOS  s  OVERLOAD  SHUT-OFF  RATIO 

(Only  for  generalized  willenborg  model) 


MS  s  MISSION  NUMBER. 


I READ  !  OPTION  FOR  INPUTTING  THE  LOAD  SPECTRUM. 

IF  IREAD  =  1,  READ  IN  SMAX,  R,  ?<  CPF 

=  2,  READ  IN  SMAX,  SMIN,  «<  CPF 
=  3,  READ  IN  SMAX,  SMEAN,  .  ft  CPF 
=  4,  READ  IN  SMEAN,  SALT,  ft  CPF 
=  5,  READ  IN  SDEL,  R,  ft  CPF 


I LAYER  s  NUMBER  OF  LOAD  LAYERS  IN  THE  MISSION. 


TITLEM  :  MISSION  IDENTIFICATION  TITLE. 


SR AT 10  :  SCALE  FACTOR  OF  THE  LOADING  SPECTRUM. 
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TABLE  2-2.  (CONTINUED) 


HOURM  :  EQUIVALENT  FLIGHT  HOURS  TO  THIS  HI SSI ON  (  IN  HOURS  ) 


SMAX  :  MAXIMUM  AXIAL  STRESS  OF  THE  LAYER.  <  UNIT  :  KSI  ) 


SMIN  :  MINIMUM  AXIAL  STRESS  OF  THE  LAYER,  i  UNIT  :  KSI  ) 


:  AXIAL  STRESS  RATIO  OF  THE  LAYER,  OR  R=SMIN/SMAX. 


CPF  :  NUMBER  OF  CYCLES  OF  THE  LAYER. 


SMEAN  :  MEAN  AXIAL  STRESS  OF  THE  LAYER.  (  UNIT  :  KSI  ) 
SMEAN  =  (SMAX+SMIN) /2. 0 


SALT  :  ALTERNATING  AXIAL  STRESS  OF  THE  LAYER.  (  UNIT  s  KSI  ) 
SALT  =  ( SMAX-SMIN ) /2.  0 


SDEL  :  AXIAL  STRESS  RANGE  OF  THE  LAYER.  <  UNIT  :  KSI  ) 
SDEL  =  SMAX-SMIN 


NPASSO  :  NUMBER  OF  BLOCKS  OF  LOAD  SPECTRUM  TO  BE  REPEATED  IF  FRACTURE 
DOES  NOT  OCCUR.  IF  NOT  SURE,  SET  NPASSO  =  0,  THEN  THE  PROGRAM 
WILL  INTERALLY  SET  NPASSO  =  1,000,000. 


NSEGT  s  TOTAL  NUMBER  OF  SEGMENTS  IN  ONE  BLOCK  OF  LOAD  SPECTRUM. 


NFLTS  :  NUMBER  OF  CONSECUTIVE  FLIGHTS  OF  SAME  MISSION  IN  ONE  SEGMENT. 


MISSION  :  MISSION  NUMBER  OF  THE  FLIGHT. 


I THRU  :  OPTION  TO  SPECIFY  FOR  THROUGH-THE-THICKNESS  OR  PART-THROUGH  CRACK. 
IF  I THRU  =0,  FOR  CORNER  CRACK, 

=  1,  FOR  THROUGH-THE-THICKNESS  CRACK. 


METHOD  :  METHOD  OF  STRESS  INTENSITY  FACTOR  CALCULATION. 

=  1,  USE  COMPOUNDING  METHOD. 

=  2,  USE  GREEN  FUNCTION  METHOD,  AND  STRESS  DISTRIBUTION 
ALONG  CRACK  GROWTH  PATH  WILL  BE  INPUT 

=-2,  USE  GREEN  FUNCTION  METHOD,  AND  AUTOMATICALLY 
CALCULATE  STRESS  DISTRIBUTION  ALONG  CRACK 
GROWTH  PATH  FROM  THE  TABLES  STORED  IN  THE  • 
PROGRAM,  WHICH  IS  A  FUNCTION  OF  OUTER-TO- 
INNER  RADIUS  RATIO  OF  LUG.  IN  THE  CASE  OF 
LUGS  WITH  BUSHINGS,  THIS  OPTION  SHOULD  BE 
EXERCISED  ONLY  WHEN  THE  RATIO  OF  LUG  OUTER 
RADIUS  TO  BUSHING  INNER  RAD I US (=P IN  RADIUS) 

IS  2.25,  SINCE  THE  SUBROUTINE  DIRCAL  CONTAINS 
DATA  ONLY  FOR  THIS  LUG-BUSHING  CONFIGURATION. 
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TABLE  2-2.  (CONTINUED) 


I BETA  :  OPTION  TO  INPUT  BETA  FACTOR'S  DIRECTLY. 

IF  I  BETA  =  O,  S  IF  WILL  BE  COMPUTED  BY  GREEN  FUNCTION  METHOD 

OR  COMPOUNDING  METHOD  DEPENDING  ON  THE  VALUE  Op 
THE  PARAMETER  '  METHOD  •' 

=  1.  DIRECT  INPUT  OF  NORMALIZED  STRESS  INTENSITY  FACTOR  * 
=-l,  AUfOMATICALLY  COMPUTE  NORMALIZED  STRESS  INTENSITY 

FACTORS  FROM  THE  TABLES  STORED  IN  THE  PROGRAM,  WHICH 
IS  A  FUNCTION  OF  OUTER-TQ- INNER  RADIUS  RATIO  OF  LUO 


I Or TYP  :  OPTION  TO  CHOOSE  THE  GREEN  FUNCTION  TYPE  <  DEFAULT  =  i  ). 

IF  I Or TYP  =  I,  USE  MODIFIED  GREEN  FUNCTION  '  FOR  LUOS  WITH  NO 
BUSHINGS  ) 

=  2,  USE  ORIGINAL  GREEN  FUNCTION  (  FOR  LUOS  WITH 
BUSHINGS  ) 


RADIUS  :  RADIUS  OF  THE  LUO  HOLE  <  UNIT  :  INCH  > 


WIDTH  :  WIDTH  OF  THE  LUG  <  =OUTER  DIAMETER  OF  LUO).  (  UNIT  :  INCH  ) 


THICK  :  THICKNESS  OF  THE  LUG.  <  UNIT  :  INCH  ) 


CRACKO  :  TMITAL  CRACK  LENGTH  'C'.  (  UNIT  :  INCH  ) 

NOTE  :  'C'  REFERS  TO  THE  CRACK  LENGTH  ON  LUG  SURFACE. 


CRACKF  s  FINAL  CRACK  LENGTH  DESIRED  IF  FRACTURE  DOES  NOT  OCCUR. 
(  UNIT  :  INCH  > 


A02CIN  s  INITIAL  FLAW  SHAPE  RATIO  OF  PART-THROUGH’ CRACK,  A/2C.  ' 

NOTE  s  'A'-  REFERS  TO  THE  CRACK  LENGTH  ALONG  LUG  HOLE  WALL. 


TITLE3  :  TITLE  FOR  METHOD  OF  COMPUTATION  OF  STRESS  INTENSITY  FACTORS. 


MBRT  :  NUMBER  OF  DIRECT  NORMALIZED  BETA  FACTORS  TO  BE  READ  IN. 


B  :  VALUE  OF  SOME  CHARACTERISTIC  LENGTH  WHICH  IS  USED  (  EXAMPLE  :  LUG 

HOLE  RADIUS  )  TO  NORMALIZE  THE  CRACK  LENGTH.  (  UNIT  :  INCH  ) 

MOTE  :  DEFAULT  =  1.0  INCH. 


COVERB  :  VALUE  OF  C/B  FOR  WHICH  THE  BETA  FACTOR  IS  BEING  PROVIDED. 


3ETAT3  :  VAuUE  OF  BETA  FACTOR  FOR  THE  ABOVE  CO VERB  VALUE. 

BETATB  =  K I / ( GROSS  SECTION  STRESS*SGRT(PIE*C> ) 


MS f G  :  NUMBER  OF  POINTS  AT  WHICH  THE  STRESS  DISTRIBUTION  IS  INPUT. 

A  NEGATIVE  VALUE  FOR  'NS 10'  MEANS  THE  SAME  AS  METHOD  =  -2- 
IM  WHICH  CASE  THE  VALUE  GF  'NSIG'  IS  DEFAULTED  TO  A  VALUE 

:.f  -  see  note  a> 
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TABLE  2-2.  (CONTINUED) 


IRSIG  :  OPT  I  !<N  TO  INCLUDE  RESIDUAL  STRESSES  I  SEE  NOTE  B>. 

IF  IRSIG  =  0.  NO  RESIDUAL  STRESSES', 

=  1,  RESIDUAL  STRESSES  ARE  ALSO  PRESENT. 

=-l,  AUTOMATICALLY  CALCULATE  THE  RESIDUAL  STRESSES 

USING  THE  INPUTS  TBR , EL , EBEL , DELD , POL  AND  POB 
(  'MS  I  O''  IS  DEFAULTED  TO  A  VALUE  OF  11  IN  THIS 
CASE  ALSO  -  SEE  NOTE  A) 


I  NR  I TEA  :  OPTION  TO  WRITE  INTERPOLATED  GREEN  FUNCTION. 

IF  IWRITEA  =  Ck  INTERPOLATED  GREEN  FUNCTION  WILL  NOT  BE  PRINTED, 
=  1-  INTERPOLATED  GREEN  FUNCTION  WILL  BE  PRINTED. 


TBR  i  BUSHING  THICKNESS  DIVIDED  BY  LOADING  PIN  RADIUS 

=  BUSHING  THICKNESS/ (LUG  INNER  RADIUS-BUSHING  THICKNE 

EL  :  YOUNG'S  MODULUS  OF  LUG  MATERIAL  (  UNIT  :  KSI  > 

EBEL  s  BUSH  I NO-TO-LUG  MATERIAL  YOUNG -'S  MODULUS  RATIO 

DELD  :  DIAMETRICAL  INTERFERENCE  LEVEL  (  UNIT  :  INCH  ) 

POL  s  PCI SSI ON'S  RATIO  OF  LUG  MATERIAL 

POB  :  POISSIONS  RATIO  OF  BUSHING  MATERIAL 

YO  :  LOCATION  AT  WHICH  THE  STRESS  IS  BEING  INPUT.  YO  =  (Y-RD/RI. 

XSIG  •  STRESS  AT  LOCATION  YO  NORMALIZED  BY  GROSS  SECTION  STRESS. 

XSIG  =  STRESS/GROSS  SECTION  STRESS. 

GROSS  SECTION  STRESS  =  P/ <2. 0*R0*TH I CKNESS* 

=  P/ ( W*TH I CKNESS  > 

RSIG  :  RESIDUAL  STRESS  VALUE  (SEE  NOTE  B).  (  UNIT  =  KSI  ) 


NOTE  A 


THE  11  POINTS  CORRESPOND  TO  DIVIDING  THE  NET-SECTION  OF  THE  LUG 
BY  10  SCiUAL  SEGMENTS  AND  THE  11  -'Y0=(  Y-RI )  /PI '  VALUES  APE 
AUTOMAT I CALL Y  CALC UL ATED . 


NOTE  0  : 


THE  RESIDUAL  STRESS  CASE  IS  APPLICABLE  TO  THROUOH-THE-TH I CKNESS 
CRACK  PROBLEMS  ONLY.  HOWEVER  FOR  CORNER  CRACK  PROBLEMS  INPUT  THE 
•MAXIMUM  STRESS  DISTRIBUTION  AND  A  STRESS  RATIO  SUCH  THAT  A 
CONSERVATIVE  LIFE  PREDICTION  COULD  STILL  BE  HADE.  (SEE  SAMPLE 
PROBLEMS  7  AMD  S) 
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SECTION  III 


SAMPLE  PROBLEMS 


In  this  section,  several  sample  inputs  and  outputs  of  the  crack  growth 
program  LUGRO  are  provided  to  aid  the  user.  Each  sample  problem  consists 
of  problem  definition,  input  and  output.  Descriptive  details  are  added  to 
the  input  and  output  to  explain  the  input  variables  and  to  interpret  the 
output.  All  the  sample  problems  illustrated  in  this  section  correspond  to 
specimens  tested  under  this  program.  Thus,  wherever  appropriate,  the  plot 
of  analytical-experimental  correlation  data  is  provided. 
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SAMPLE  PROBLEM  #1 


Material 

Aluminum 

VRi 

2.25 

R. 

1 

0.75 

Thickness 

0.50  Inch 

*0 

6  Ksi 

R 

0.1 

Crack  Type 

Through-the-Thickness 

c 

0.025  Inch 

Crack  growth  analysis  for  the  above  lug  and  crack  geometries  and  loading 
condition  is  conducted  in  this  sample  problem  using  Forman’s  crack  growth 
rate  equation  and  the  following  given  ft  factors. 


No. 

c/(Rn  -  R. )* 

P 

1 

0.0 

5.957 

2 

0.1 

3.899 

3 

0.2 

3.126 

A 

0.3 

2.804 

5 

0.  A 

2.500 

6 

0.5 

2.422 

7 

0.6 

2.364 

8 

0.7 

2.463 

9 

0.8 

2.646 

10 

0.9 

3.308 

*Rq  -  R  =  B  =  0.9375  Inch 


Analytical-experimental  correlation  result  is  given  in  Figure  3-1. 
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Line  1:  Title 

Line  2:  IEQU  =  2,  Select  Forman  Equation 


Lines  3, 4, 5, 6:  Material  Property 

data  AK  vs.  da/dN 

AK.  .  .  K  and  F 
threshold  c  ty 


-0.  1 

74 _  o 


Line  7:  IRETAR  =  0,  no  retardation 


100.0 

Line  16:  ITHRU  =  1,  IBETA  =  1 

Through- the-Thickness  cracky 
Direct  input  of  8  factors 


Lines  9,10,11, 
12,  13: 
Constant 

Amplitude  Loadin 
Spectrum  Definiti 


O. o  |  Lines  17,  18:  Lug  and  Crack  Geometries 


Lines  20,21,  22:  Input  of  8  factors 


25 


0 :*  <r  M  *N 
C<  If'  N  CO  CO 
01  04  04  04  W 


N  O’-  NOc 
—  ©  0  -6  — » 
^  o  cc*  <n  n 
©  0/  <•  in  IT* 

CO  O  04  CO  rs* 
—  ^  (4  <SJ  00 


0-  -  V>  O'*  0 

<i  O  «  ■?  « 
C*  O  00  ©  <• 

c  o>  r*  >s 

(V*  CO  .v.  «T  nO 


N 

04 


CO 

CO 

u. 

© 


©  ©  “  o  c 

UJ 

cc 

© 

-J 

©  z  ©  ©  c 

c 

u. 

LV  ONO-O 

coco 
u  cc 

U'O 
>  X 
o 


©  o 

©  ~ 

Id 

o© 

w  © 
© 
© 


CO 

GO 
H-  H 

H-  K 

© 

z  z 

-1 

Ui  Ui 

**ON{KO 

Ui 

_i 

<r  <r  — «  — 

<L  <Z 

cKci’O 

> 

>  > 

ro  c  -0  <r  o'> 

M  H-< 

C  ©  O  —  r- 

z 

”<  “l 

o 

c  c- 

*-< 

UJ  u. 

H- 

U* 

CO  if* 

Hi 

CO 

CO  © 

</>  o> 

Ilf 

©  c 

fi  «r  C  <T  N 
©  U>  -  If*  CO 

z 

Cl  l. 

»r  ©  <>  —  ro 

c 

U  Ui 

H 

z  z 

0.»  04  oi  Q- 

«N  10  ^  0«  © 

Ui 

©  o 

«N  CJ  0-i  0,  C4 

© 

•1 

M 

u 

»- 

c* 


SAMPLE  PROBLEM  it  2 


This  sample  problem  is  identical  to  the  previous  sample  problem, 
except  that  the  /S  factors  are  obtained  from  an  already-available  table  in 
the  program,  instead  of  providing  them  as  input.  This  is  accomplished  by 
specifying  the  value  of  'IBETA’  as  -1  in  the  input. 

Since  the  factors  input  in  the  sample  problem  //I  are  essentially  the 
same  as  those  stored  in  the  table,  the  results  of  crack  growth  history  are 
also  the  same.  The  slight  difference  in  &he  solutions  (<.15%)  is  due  to 
computational  round-offs. 


1  :  ABPLC46  AND  ABPLC93 


!d 

0 

0 

l 
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4 
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0 

145 
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5 . 0 

0 
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900 
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000 
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30.0 

353 
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14 
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15 

45.0 

6520 

000 

16 

60 . 9 
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17 
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0 . 9 

74. 

9 

IS 
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19 
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1 

20 

CONSTANT 
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DE  WITH 

SIGMA (MAX )=  6  KS 

H 

2.  J. 

1.0 

1 . 0 

22 

6 . 0 

0.  1 

100. 

o 

24 

w 

1 

25 

1 

1 

Q> 

do 

i 

± 

0 

_  Line  16: 

’Z*  T 

0. 

7500 

Tj 

375  C 

.  5000 

‘JW 

0. 

0250 

m 

375 

I  AND  R=0. 1 


IBETA  =  -1,  Calculate  3 
factors  from  the  tables 
stored  in  the  program 
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SAMPLE  PROBLEM  03 


This  sample  problem  also  analyzes  the  same  lug  and  crack  geometries 
and  loading  condition  as  the  two  previous  sample  problems.  However,  in 
this  case,  the  Green's  function  method  is  used  for  computation  of  stress 
intensity  factors.  The  following  stress  distribution  (normalized  by  the 
gross-section  stress,  oq )  is  input  as  a  function  of  (y  -  R^/R^ 


No. 

Cy  -  RiJ/Rj 

a/aQ 

1 

0.0000 

5.3190 

2 

0.0625 

4.2390 

3 

0.1875 

3.0938 

4 

0.3125 

2.3490 

5 

0.4375 

1.9395 

6 

0.5625 

1.6268 

7 

0.6875 

1.3950 

8 

0.8125 

1.1700 

9 

0.9375 

0.9720 

10 

1.0625 

0.7200 

11 

1.1875 

0.4073 

No  input  is  provided  for  the  parameter  'IGFTYP*.  The  program  thus 
defaults  the  value  of  'IGFTYP'  to  1,  selecting  the  modified  Green's  func¬ 
tion  to  calculate  the  stress  intensity  factors.  Comparison  of  results  of 
this  sample  problem  with  the  previous  two  sample  problems  shows  a  crack 
growth  life  difference  of  about  2%,  which  arises  from  slight  differences  in 
the  stress  intensity  factors  in  the  two  cases. 


ABPLC46  AND  ABPLC93 
2  0  0  1 
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°  O 
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0 
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1 . 9395 
1 . 6263 
1 . 3950 
1. 1700 
0.9720 
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0. 4073 


100.0 

/ —  Line  16: 


0. 5000 


METHOD  =  2,  IGFTYP  =  blank. 
Use  Green’s  function  method, 
select  modified  Green's 
function. 


Lines  23,  25:  Stress  distribution 
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PREDICTED  CRACK  GROWTH  HISTORY  O.INIFOPM  THR'I-CRACK  REGION) 


SAMPLE  PROBLEM  #4 


Instead  of  inputting  the  normalized  stress  distribution  as  in  the 
previous  sample  problem,  it  can  also  be  obtained  from  an  already-available 
table  in  the  program.  This  can  be  accomplished  by  specifying  'METHOD'  =  -2 
or  using  a  negative  value  for  'NSIG'.  This  is  illustrated  in  the  two  sets 
of  input  data  provided  for  this  sample  problem.  The  output  data  for  the 
two  sets  of  input  data  will  be  identical. 
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SAMPLE  PROELEM  #5 


Material 

R  /R. 
o  1 

R. 

l 

Thickness 


Crack  Type 

.. 


Steel 

3.0 

0.75  Inch 
(J.5  Inch 
14  Ksi 
0.1 

Corner 
0.025  Inch 
0.5 


This  sample  problem  analyzes  and  predicts  the  growth  behavior  of  a 
corner  crack  using  the  Green’s  function  method.  ’ITHRU’  =  0  specifies  that 
this  is  a  corner  crack  problem.  A  value  of  0.5  is  input  for  ’A02CIN*  to 
specify  that  the  shape  of  the  initial  corner  crack  is  quarter-circular. 
The  output  contains  three  regions  of  crack  growth,  namely  corner  crack 
growth  until  the  crack  breaks  through  the  thickness,  transitional  crack 
growth  to  through-the-thickness  crack  shape,  and  subsequent  through-the- 
thickness  crack  growth  until  failure. 

Analytical-experimental  correlation  result  for  this  sample  problem  is 
given  in  Figure  3-2. 
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PREDICTED  PART-THROUGH  CORNER  CRACK  GROWTH  HISTORY 
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NOTE  t  ONE  PASS  IS  EQUIVALENT  TO  100.0000  CYCLES 

ONE  PASS  IS  EQUIVALENT  TO  1.0000  FLIGHT  HOURS 


SAMPLE  PROBLEM  #6 


This  sample  problem  illustrates  the  use  of  compounding  method  to 
analyze  a  crack  problem  by  specifying  'METHOD’  =  1.  The  lug  and  crack  geo¬ 
metries  and  loading  condition  are  the  same  as  the  sample  problem  #5.  In 
the  output,  only  one  crack  length  (the  front  surface  crack  length  c)  is 
printed.  Since  a/c  ratio  is  assumed  as  constant  in  the  compounding  method, 
the  crack  length  along  the  lug  hole  wall  or  the  back  surface  crack  length 
can  be  calculated  using  the  front  surface  crack  length,  a/c  ratio  and  the 
thickness  of  the  lug. 

The  computed  solution  using  the  compounding  method  is  also  included  in 
Figure  3-2  (sample  problem  # 5). 
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SAMPLE  PROBLEM  #7 


Material 

Aluminum 

VRi 

2.25 

Ri 

0.75  Inch 

Thickness 

0.5  Inch 

Crack  Type 

Through-the-Thic  knes  s 

co 

0.025  Inch 

Loading 

As  Described  Below 

This  sample  problem  illustrates  the  analysis  of  a  through-the-thick- 

ness  crack  problem  with  residual  stress.  In  this  case,  the  lug  is  loaded 

with  a  far-field  stress  ( a  )  of  15  ksi  and  a  stress  ratio  (R)  of  0.3.  For 

o 

this  loading  condition,  the  uncracked  lug  undergoes  plastic  yielding.  An 
elasto-plastic  stress  analysis  was  conducted  tc  obtain  the  stress  distribu¬ 
tion  at  the  maximum  and  minimum  (corresponding  to  R  =  0.5)  stress  levels. 
The  resulting  stress  distributions  are  presented  in  Figure  3-3.  The  stress 
distribution  at  the  minimum  stress  level  is  considered  as  the  residual 
stress  and  the  difference  in  the  stress  distributions  between  maximum  and 
minimum  stress  levels  is  considered  as  the  alterating  stress.  This  stress 
distribution,  without  normalizing  by  the  far-field  stress  (o^),  is  speci¬ 
fied  as  shown  in  the  input  data.  Since  the  stress  distribution  is  not 
normalized,  values  of  1.0  and  0.0  are  input  for  'SMAX*  and  'R'  in  the 
spectrum  layer  input. 

Comparison  of  analytical  and  experimental  results  is  presented  in 
Figure  3-1*.  The  present  output  data  correspond  to  the  curve  labelled 
PLASTIC  STRESS. 
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Line  12:  SMAX  =  1.0,  R  =  0.0 


Line  16:  ITHRU  =  1,  through-the- 
thickness  crack 

Line  23:  IRSIG  =  1,  residual  stress 
present 


Line  25*:  Stress  distribution 
(cyclic  and  residual 
stresses) 
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CRACK  LENGTH,  c,  INCH 


SAMPLE  PROBLEM  #8 


In  the  input  instruction  (Section  II),  Note  B  states  that  the  residual 
stress  case  is  applicaole  to  the  through-the-thickness  crack  problem  only. 
Sample  problem  #7  illustrated  the  analysis  of  a  through-the-thickness  crack 
problem  with  residual  stress.  Note  B  also  states  that,  for  analyzing  cor¬ 
ner  crack  problems  with  residual  stress,  input  the  maximum  stress  distribu¬ 
tion  and  a  stress  ratio  such  that  a  conservative  life  prediction  can  still 
be  made.  Using  the  data  of  sample  problem  #7*  the  maximum  stress  distribu¬ 
tion  and  a  conservative  R  are  calculated  as  follows: 


(y-R^/R 

i  ^alt 

°Yes 

°max  ^alt. 

+  a  ) 
res 

R  _  * co 

eff  "  ^max 

(KSI) 

(KSI) 

(KSI) 

0.0000 

34,70 

32.10 

66.80 

0.48054 

0.0625 

30.97 

28.22 

59.19 

0.47677 

0.1875 

22.84 

22.54 

45.38 

0.49669 

0.31?5 

17.70 

18.26 

35.96 

0.50779 

0.4375 

14.79 

15.40 

30.19 

0.51010 

0.5625 

12.46 

12.86 

25.32 

0.50790 

0.6875 

10.73 

11.06 

21.79 

0.50757 

0.8125 

9.02 

9.30 

18.32 

0.50764 

0.9375 

7.53 

7.80 

15.33 

0.50881 

1.0625 

5.66 

5.95 

11.61 

0.51249 

1.1875 

3.31 

3.61 

6.92 

G. 52168 

A  conservative  R 

will  be  the  minimum  value  of  R 

in  the  above  table. 

The  maximum  stress  distribution 

and  the  conservative  R 

are  then  specified 

as  shown 

in  the  input 

data. 

Also 

note  in  the  output  < 

data  of  sample 

problems 

1  through  6,  the 

failures 

are  due  to  net-section 

yielding.  In 

sample  problems  7  and  8,  the 

failures 

are  due  to 

the  exceedence  of  critical  stress  intensity  factor 

(fracture  toughness)  values. 

Predicted  and  experimental  data  are  compared  in  Figure  3-5.  The  curve 
labelled  PLASTIC  STRESS  corresponds  to  present  output  data. 
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Line  12:  SMAX  =  1.0,  R  =  0.47677 


Line  16:  ITHRU  =  0,  corner  crack 
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NOTE  >  ONE  PASS  IS  EQUIVALENT  TO  100.0000  CYCLES 

ONE  PASS  IS  EOUI VALENT  TO  1.0000  PLIGHT  HOURS 


Figure  3-5.  Corner  Crack  Growth  Data  and  Prediction, 

Aluminum  Lug,  Ro/R.=2.25,  <To=15  KSI,  R=0.5 


SAMPLE  PROBLEM  #9 


In  all  the  sample  problems  presented  thus  far,  only  simple  constant 
amplitude  loading  i3  illustrated.  This  sample  shows  the  input  for  a 
slightly  complex  block  spectrum  loading.  The  schematic  ar.d  details  of 
block  spectrum  loading  are  given  in  Figure  3-6  and  Table  3-1.  The  method 
of  inputting  the  block  spectrum  is  shown  in  the  input  data. 

The  above  block  spectrum  will  introduce  spectrum  retardation  effects. 
The  Generalized  Willenborg  model  CIRETAR1  =  1  and  'IMODEL*  =  2)  is  used  in 
this  sample  problem  to  account  for  the  crack  growth  retardation  effect. 
Values  of  0.0  and  3.5  are  used  in  the  input  for  the  parameters  'XKMXTH1 
(maximum  threshold  stress  intensity  factor)  and  'SOS*  (overload  shut  off 
ratio),  respectively.  These  values  correspond  to  a  value  of  0.4  for  the 
factor  0  of  the  Generalized  Willenborg  equation.  The  plastic  zone  is  cal¬ 
culated  assuming  a  plane  stress  condition  CIPLANE'  =  1). 

Figure  3-7  shows  the  analytical-experimental  data  correlation  for  this 
sample  problem.  The  labels  mean  the  following. 

[A]  -  Hsu  retardation  model  is  used 

[B]  -  Generalized  Willenborg  model  is  used  (Present  output  data) 

[C]  -  Willenborg  model  is  used 

[D]  -  No  retardation  effect  is  considered. 
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STRESS  RATIO,  R  =  0. 1 


NUMBER  OF  CYCLES 

Figure  3-6.  Schematic  of  Block  Spectrum  Loading. 


TABLE  3.1.  DETAILS  OF  BLOCK  SPECTRUM  LOADING 


STRESS  DATA 


LOAD  FACTOR 


SYMBOL 


MAX.  STRESS 


CYCLE 

DATA 

SYMBOL 

NO.  OF 
CYCLES 

Hi 

950 

N2 

650 

-3 

450 

«« 

250 

136 

44 

*>7 

15 

B8 
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Nft 
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UNIT  LOAD  FACTOR  =  6  KSI;  R 
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ABPLS/6 


4: 

4 . 0 

0.  145 

5: 

5.0 

0.573 

6: 

6.  0 

2 . 050 
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i  • 
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1 0 . 0 
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1 3 . 200 
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23.300 

1 1 : 

20 . 0 

53. 900 
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25 . 0 

x  34 . 000 

IS: 

30 . 0 

353. 000 

14: 

40.0 

2850. 000 

15: 

45.  0 

6520 . 000 

16: 

60.  94 

0.  1 

17: 

0 . 0 

60.  94 

-0.  1 

74.9 


BLOCK  Si 


PSCTRl 
1 . 0 
2.7 


UM  DEFINITION 


1.0 

0.  1 

950 . 0 

0.  1 

650.0 

0.  1 

450.0 

0  1 

250 . 0 

v  •  i. 

0.  1 

136.0 

0.  1 

44.0 

0, 1 

15.0 

0.  1 

4.0 

0.  1 

1.0 

Line  7: 


IRETAR  =  1,  IMODEL  =  2, 
Account  retardation  effect 
with  Generalized  Wilienborg 
model 


Line  8:  IPLANE  =  1, 

XKMXTH  =0.0, 
SOS  =  3.5 


Line  12:  Block  Spectrum  Loading 
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CRACK  LENGTH.  C/  INCH 


PASSES 


Figure  3  7.  Through-the-Thickness  Crack  Growth  Data  and  Prediction,  Aluminum 

iuS»  Rd/i=2'2~’’  E^ock  Spectrum  Loading,  <j  =7.5  KSI 

omax 


SAMPLE  PROBLEM  #10 


li- 

a 
-** 
m  ; 


This  example  corresponds  to  prediction  of  corner  crack  growth  behavior 
in  a  lug  subjected  to  a  much  more  complex  spectrum  of  flight-by-flight 
leading.  This  spectrum  is  typical  for  a  cargo  aircraft  (C-5).  This  spec¬ 
trum  inclines  thirteen  different  missions  (0,  1  through  12).  The  details 
of  the  missions  are  provided  in  Table  3-2.  In  this  table,  a  special  load¬ 
ing  cycle  is  defined  in  terms  of  N/FLT  =  0.1,  for  example  in  Mission  1. 
This  simply  means  that  the  particular  load  is  applied  once  in  ten 
occurrences  of  this  mission.  When  Mission  1  includes  this  load  it  is 
referred  to  as  Mission  1*.  In  the  input  for  this  problem,  Mission  1  and 
Mission  1*  are  treated  as  separate  missions.  The  sequence  of  missions  in 
one  pass  (mission  mix)  of  the  spectrum  is  given  in  Table  3-3.  One  pass 
consists  of  120  missions  and  one  mission  has  the  same  meaning  as  a  flight. 

In  this  sample  problem,  the  Hsu  retardation  model  CIRETAR'  =  1  and 
'IMODEL'  =3)  is  used  and  the  yield  zone  is  calculated  using  the  plane 
strain  condition  ( ' IPLANE'  =  2). 

Analytical  and  experimental  results  are  compared  in  Figure  3-8. 
Labels  [A],  [B],  [C]  and  [D]  mean  the  same  models  as  given  in  sample 
problem  #9. 
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TABLE  3-2 •  MISSIONS  DEFINITION  FOR  CARGO  SPECTRUM 


MAX 

MIN 

STRESS 

STRESS 

K3I 

KSI 

M/PLT 

MS  0 

19.108 

co 

T  — 

• 

1 

1  .00 

MIS  1 

18.083 

-1.808 

1  .00 

4.423 

2.423 

1998. 00 

13.492 

9.492 

193.00 

15. 106 

9.106 

24.00 

15.863 

7.363 

4,00 

16.909 

5.572 

1 .00 

18.083 

2.602 

.10 

MIS  2 

17.987 

-1.799 

1  .00 

5.080 

3.080 

2273. CO 

13.324 

9.824 

204.00 

15.189 

9.189 

25.00 

15.908 

7.908 

5.00 

16.899 

5.494 

1  .00 

17.987 

2.384 

.10 

MIS  3 

16.751 

-1.675 

1 .00 

4.888 

2.888 

2891 .00 

14.139 

10.139 

227.00 

15.304 

9-304 

25.00 

15.918 

7.918 

4.00 

16.751 

5.302 

1 .00 

IMS  4 

16.695 

-i . 669 

1  .00 

6.157 

4.157 

2827.00 

14.755 

10.755 

209.00 

15.562 

9.562 

23.00 

16.065 

8.065 

4.0C 

16.695 

5.536 

1  .00 

MS  5 

17.952 

-1.795 

1  .00 

4.737 

2.737 

2099.00 

14.329 

10.329 

175.00 

15.420 

9-420 

21  .00 

16.101 

8.101 

4. CO 

16.936 

5.964 

1  .00 

17.952 

2.388 

.10 

*  N/FLT  =  0,1 -Means  the  application  of 

this  load  once  in 

~er_  occurrence: 

of  thi3  Ms sion  (refer  to  sequence  of 

Mssions  table). 
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TABLE  3-2. 

MISSIONS 

DEFINITION  FOR  CARGO  SPECTRUM 

(CONTINUED) 

MAX 

STRESS 

KSI 

MIN 

STRESS 

KSI 

N/FLT 

MIS  6 

16.778 

-1 .678 

1  .00 

6.313 

4.313 

5082.00 

15.372 

11.372 

204.00 

15.904 

9-904 

21.00 

16.521 

8.321 

3. CO 

16.778 

6.050 

1  .00 

MIS  7 

17.861 

-1.786 

1.00 

5.628 

3.628 

2100.00 

111.  634 

10.634 

164.00 

15.592 

9.592 

19.00 

16.232 

8.232 

3.00 

16.916 

6.358 

1 .00 

17.861 

3-335 

.10 

MIS  S 

17.518 

-1 .752 

1 .00 

6.4l4 

4,4l4 

2714.00 

15.401 

n.401 

187.00 

15.941 

9.941 

19.00 

16.373 

3.373 

3.00 

16.820 

6.253 

1.00 

17.518 

2.991 

.10 

MIS  9 

17.761 

-1 .776 

1 .00 

6.5H4 

4.544 

2260.00 

15.197 

11.197 

161 .00 

15.853 

9.858 

17.00 

16.406 

8.406 

2. CO 

16.913 

6.694 

1.00 

17.761 

3.686 

.10 

MIS  10 

TABLE  3-2.  MISSIONS  DEFINITION  FOR  CARGO  SPECTRUM  (CONTINUED) 


MAX 

MIN 

STRESS 

STRESS 

ESI 

ESI 

i'l/FLT 

MIS  11 


14.546 

-1 .808 

1 .00 

15.523 

.000 

1  .00 

15.523 

.000 

1 .00 

<5.523 

.000 

1  .00 

18.083 

.000 

1 .00 

15.523 

.coo 

1.00 

15.523 

.000 

1.00 

15.523 

.000 

1.00 

16.191 

.000 

1 .00 

16.191 

.000 

1 .00 

3.774 

1 .774 

5418.00 

12.842 

8.842 

419.00 

1 4 . 546 

8.546 

51.00 

15.523 

7.523 

9.00 

16.191 

6.191 

2.00 

16.913 

4.425 

1.00 

18.083 

1 .151 

.10 

15.650 

-1 .849 

1 .00 

15.650 

.000 

4 .00 

15.650 

.000 

1.00 

15.650 

.000 

1  .GO 

18.487 

-1 .849 

1 .00 

15.650 

.000 

1  .00 

16.323 

.000 

1 .00 

16.323 

.000 

1 .00 

16.323 

.000 

1 .00 

3.848 

1.848 

6512. CO 

13.112 

9.112 

556.00 

14.771 

3.771 

69.00 

15.650 

7.650 

12.00 

16.323 

6.323 

5.0  0 

17.267 

4.096 

1  .00 

18.487 

.919 

.10 
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TABLE  3-3.  ONE  PASS  OF  SEQUENCE  OF  MISSIONS  OF  CARGO  SPECTRUM 


SEQUENCE 

MISSION 

SEQUENCE 

MISSION 

SEQUENCE 

MISSIOJ' 

NO. 

NO. 

NO. 

NO. 

NO. 

NO. 

1 

7 

41 

2 

81 

7 

2 

8 

42 

5 

82 

12 

3 

1 

43 

12 

83 

11 

4 

2 

44 

7 

84 

9 

5 

12* 

45 

2* 

85 

7 

6 

7 

46 

3 

86 

3* 

7 

5 

47 

1 

87 

2 

3 

11 

43 

9 

88 

5 

9 

8 

49 

4 

89 

1 

10 

9 

50 

7 

90 

12 

11 

*7 

51 

12 

91 

7 

12 

1 

52 

3 

92 

10 

13 

12 

53 

11 

93 

8 

14 

2 

54 

5 

94 

11 

15 

7* 

55 

7* 

95 

3 

16 

3 

56 

1 

96 

7* 

17 

11 

57 

2 

97 

2 

13 

5 

53 

12 

98 

1 

19 

1 

59 

8 

99 

12 

20 

7 

60 

7  ■ 

100 

5 

21 

12 

61 

0 

101 

7 

22 

4 

62 

3 

102 

3 

23 

2 

63 

5 

103 

9 

24 

7 

64 

1 

104 

2 

25 

5* 

65 

7 

105 

7 

26 

6 

66 

12* 

106 

11* 

27 

11 

67 

3 

107 

12 

23 

1 

68 

2 

108 

1 

29 

12 

69 

11 

109 

8 

30 

7 

70 

9 

110 

5 

31 

3 

71 

7 

111 

7 

32 

2 

72 

1 

112 

2 

33 

5 

73 

12 

113 

12 

34 

7 

74 

3 

114 

3 

35 

1* 

75 

7 

115 

7 

36 

9 

76 

5* 

116 

9* 

37 

12 

77 

2 

117 

1 

38 

3 

78 

4 

118 

11 

39 

11 

79 

1 

119 

4 

40 

7 

80 

3 

120 

12 

"Missions  with  application  of  once  in  ten  occurrances  loads  (i.e..  Loads 
with  N/FLT  =0.1) 
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Line  7:  IRETAR  =  1,  IMODEL  = 
Account  retardation 
effect  with  Hsu  model 


g 


i 

SBPLS56  AND 

SBPLS54 

2 

2  0 

0  0 

**» 

13  0 

4 

4.0 

0.0537 

5 

5.  0 

0.  106 

6 

7.0 

0.247 

7 

9.0 

0.  482 

o 

11.0 

0.833 

*;/ 

13.0 

1.310 

10 

16.0 

2.  250 

11 

20.0 

3.  930 

12 

24.  0 

6.010 

13 

28.0 

8.350 

14 

35.0 

12.800 

15 

50.0 

29. 000 

16 

58.0 

57 . 700 

17 

224.7 

0.5 

-0.  1 

18 

0.0 

224.7 

179.7 

19 

20 

21 

_ 

. 

Y  2~ 

T”  ' 

—  -  - 

22 

MISSION  0 

23 

1.0 

1.0 

24 

19. 108 

-1.184 

1.0 

25 

2  2 

6 

26 

MISSION  1 

27 

1.0 

1.0 

28 

18.083 

-1.308 

1.0 

29 

4.423 

2.  423 

1998.0 

30 

13.492 

9.492 

193.0 

31 

15. 106 

9.  106 

24.0 

82 

15.863 

7.863 

4.0 

16.909 

5.572 

1.0 

34 

7 

35 

MISSION  1* 

36 

1 . 0 

i.O 

37 

18. 083 

- 1 . 808 

1.0 

OjZ; 

4.423 

2.423 

1998.0 

39 

13.492 

9.492  ‘ 

193.0 

40 

15. 106 

9.  106 

24.0 

41 

15.363 

7.863 

4.0 

42 

16.909 

5.572 

1.0 

43 

18.083 

2.  602 

1.0 

44 

4  2 

45 

MISSION  2 

46 

1.0 

1.0 

47 

17.987 

-1.799 

1.0 

43 

5 . 080 

3.080 

2273.0 

49 

13.824 

9.  324 

204 . 0 

50 

15. 189 

9.  189 

25 . 0 

Line  8:  IPLANE  =  2 


Lines  9,10,11,  12*,  13: 
Flight-by-flight 
cargo  spectrum 
definition  beginning 


m 


51: 

15.908 

7.  903 

5.0 

52: 

16.899 

5.  494 

1 . 0 

53: 

5  2 

7  ■ 

54: 

55: 

MISSION  2* 
1.0 

1 . 0 

56: 

17.987 

-1.799 

1.0 

57: 

5.080 

3.080 

2273.0 

58: 

13.824 

9.824 

204.0 

5?: 

15. 189 

9.  189 

25.0 

60: 

15.908 

7.908 

5.0 

61 : 

16.899 

5.  494 

1 . 0 

62: 

17.987 

2.384 

1.0 

63: 

/z*  2 

6 

64: 

65: 

MISSION  3 
1.0 

1.0 

66: 

16.751 

-1.675 

1.0 

67: 

4 . 888 

2.888 

2891.0 

68: 

14. 139 

10.139 

227.0 

6'?: 

15.304 

9. 304 

25.0 

70: 

15.918 

7.918 

4.0 

71: 

16.751 

5.302 

1.0 

72: 

7  2 

6 

78* 

74: 

MISSION  4 
1.0 

1.0 

75: 

16.695 

-1.669 

1.0 

76: 

6.  157 

4. 157 

2827.0 

77: 

14.755 

10.755 

209.0 

78: 

15.562 

9.562 

23.0 

79: 

16.065 

8.065 

4.0 

80: 

16.695 

5.536 

1.0 

81* 

8  2 

6 

OO  • 

—  •*—  ■ 

83: 

MISSION  5 
1.0 

1.0 

84: 

17.952 

-1.795 

1 . 0 

•:<c  • 

■j  • 

4.737 

2.737 

2099.0 

86: 

14.329 

10.329 

175.0 

87: 

15.420 

9.420 

21.0 

88: 

16. 101 

8.  101 

4.0 

89: 

16.936 

5.964 

1.0 

90: 

9  2 

7 

91: 

92: 

MISSION  5* 
1.0 

1.0 

93: 

17.952 

-1.795 

1.0 

94: 

4,737 

2.737 

2099.0 

95: 

14.329 

10.329 

175.0 

96: 

15.420 

9 . 420 

21.0 

97: 

16. 101 

8.  101 

4.0 

OO  • 

16.936 

.5.964 

1.0 

go  • 

.  s  * 

17.952 

2.  888 

1.0 

1 00 : 

10  2 

92 


C3EC3TS  WH  JL (UM.kiSJL'H  L’.V iTVA *.%.  1.--  i-V  JJV.  I.KIL SJIi uU «  .Tv 


:z,l  >  *  J  ^  J, yvv ,-jlv.M- , _* AT I’VJ^-'LW  Wl ; 


102  - 

t.  0 

1.0 

103: 

16. 773 

--1 . 673 

1.0 

104: 

6.313 

4.313 

3032. 0 

1 05 : 

15.372 

1  i . 372 

204 .  (j 

106: 

15.904 

.  904 

21.0 

1 07 : 

16.321 

3.321 

3 . 0 

103: 

16.773 

6.  050 

1 . 0 

109: 

11  2 

6 

1 10: 

MISSION  7 

111: 

1.0 

1 . 0 

?  ■*.  •“  ; 

17.361 

-1.736 

1 . 0 

1 13: 

5.623 

*!  a  <£*2‘Z' 

2100. 0 

114: 

14.634 

10.634 

164.0 

115: 

15.592 

9.592 

1 9 . 0 

1 16: 

16.232 

y .  v^2 

3 .  0 

117: 

16. 916 

6.353 

1 . 0 

1  IS: 

12  2 

7 

119: 

MISSION  7* 

120: 

1 . 0 

1.0 

121: 

17.861 

- 1 . 736 

1 . 0 

122: 

5 . 623 

3.623 

2100.0 

123: 

14.634 

10.634 

16*1.0 

124: 

15.592 

592 

19.0 

125: 

16.232 

— *  • 

3 .  0 

126: 

16.916 

6.  353 

1.0 

127: 

17.361 

3.  335 

1.0 

i  23 : 

13  2 

6 

129: 

MISSION  3 

130: 

1 . 0 

1.0 

131 : 

17.513 

-1.752 

1.0 

132: 

6.414 

4.414 

2714.0 

133: 

15.401 

11.401 

137.  0 

134: 

15.941 

9.941 

1 9 . 0 

135: 

16. 373 

3.373 

3.  0 

136: 

16. 320 

OS": 

1 . 0 

1 37 : 

14  2 

7 
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MISSION  3* 

139: 
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1 . 0 

140: 

17.518 

-1.752 

1.0 

141 : 

6.414 
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2714.0 

142: 

15. 401 

11.401 

137.0 

143: 

15. 941 

9.  941 

1 9 . 0 

144: 

16. 373 

3.373 

3.  0 

145: 

16. 320 

6-253 

1 . 0 

146: 

17.513 

2.  991 

1 . 0 

147: 

15  2 

6 

143: 

MISSION  9 

149: 

i  .  0 

1 . 0 

1 50 : 

17.761 

-1.776 

1 . 0 
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rrtitURE  rjur.  to  met  section  yieU'1  (  ■=-tre^>  -  2^.7 *  k-m  > 


SAMPL"  PROBLEM  #11 


Analysis  of  a  lug  with  an  interference-fit  bushing  is  illustrated  in 
this  sample  problem.  Details  of  the  material  and  geometric  parameters  are: 


Lug  Material 
Bushing  Material 
Lug  Outer  Radius,  Rq 
Lug  Inner  Radius,  R^ 

Lug  Thickness 
Bushing  Thickness,  t 

3 

Young's  Modulus  of  Lug  Material,  E^ 
Young's  Modulus  of  Bushing  Material,  E 

O 

Poisson's  Ratio  of  Lug  Material,  V L 
Poisson's  Ratio  of  Bushing  Material,  V 

D 

Loading  Pin  Radius  (=  Lug  Inner  Radius- 
Bushing  Thickness),  r^ 

Vi 

Vl 

Diametrical  Interference,  6 „ 


Crack  Type 


Steel 

Steel 

1.6875  Inch 
0.84  Inch 
0.5  Inch 
0.09  Inch 
3  x  10**  Ksi 
3  x  104  Ksi 


0.75  Inch 


0.008  Inch 
14  Ksi 


Through-the-Thickness 
0.025  Inch 


Four  sets  of  input  data  are  presented  for  this  problem.  In  the  first 
data  set,  the  residual  stress  and  the  normalized  stress  distribution  across 
the  net  section  are  calculated  from  the  already-available  tables  in  the 
program  by  specifying  'METHOD'  =  -2  and  'IRSIG'  =  -1.  In  the  second  data 
set,  the  residual  stress  and  the  normalized  stress  distribution  (obtained 
using  the  concentric  cylinder  equations  and  the  finite  element  method)  are 
input  directly.  In  the  third  data  set,  the  residual  stress  is  calculated 
automatically  ( 'IRSIG'  =  -1)  and  the  normalized  stress  distribution  is 
input  directly.  In  the  fourth  data  set,  the  residual  stress  is  input 
directly  and  the  normalized  stress  distribution  is  obtained  from  the 


already-available  tables  ('METHOD'  =  -2).  In  the  fourth  data  set,  note 
that  only  tn/r.  CTBR')  and  ED/E.  ('EBEL')  are  needed  to  interpolate  the 

D  1  D  L 

dara  from  the  tables.  Values  for  parameters  'EL',  'DELD',  'POL'  and  'POB' 
are  not  needed  for  such  an  interpolation,  and  will  be  ignored  by  the 
program  if  specified,  without  influencing  the  output  results. 

The  above  four  input  data  sets  show  the  flexibility  with  which  the 
stress  data  can  be  input.  All  the  four  input  data  sets  are  equivalent  and 
thus  only  one  set  of  output  data  is  included  for  this  sample  problem.  The 
other  three  sets  of  output  data  would  be  the  same,  except  for  slight 
differences  in  the  results  due  to  computational  round-offs. 

The  far-field  loading  is  applied  with  a  stress  ratio  of  0.1.  However, 
the  operating  stress  ratio  across  the  net  section  of  the  lug  will  be  dif¬ 
ferent  due  to  the  presence  of  residual  stress.  In  this  case,  the  operating 
stress  ratio  ranges  from  0.32  at  the  inner  radius  location  to  0.64  at  the 
outer  radius  location.  Thus,  for  this  problem,  crack  growth  rate  data 
corresponding  to  a  stress  ratio  0.5  are  input,  although  the  applied  far- 
field  stress  ratio  is  0.1. 

A  value  of  2  has  been  input  for  parameter  'IGFTYP' ,  selecting  the 
original  (unmodified)  Green's  function  for  the  computation  of  stress 
intensity  factors  in  this  problem.  The  use  of  the  original  (unmodified) 
Green's  function  ('IGFTYP'  =  2)  is  recommended  for  lugs  with  interference- 
fit  bushings.  The  modified  Green’s  function  ('IGFTYP'  =  1,  the  default 
value)  should  be  used  for  lugs  with  no  bushing,  since  the  modification  of 
Green's  function  accounts  for  the  pin  bearing  pressure  distribution 
variation  with  respect  to  crack  length. 

Comparison  of  experimental  and  analytical  data  of  this  sample  problem 
is  given  in  Figure  3-9. 
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NOTE  t  ONE  F  ASS  IS  E OU1 VALENT  TO  100.0000  CYCLES 

ONE  PASS  IS  EQUIVALENT  TO  1.0000  FLIGHT  HOURS 


CRACK  LENGTH,  c,  INCH 


Figure  3-' 


.  Through-the-Thickness  Crack  Growth  Data  and  Prediction,  Steel  Lug 
with  Steel  Bushing,  R  /r.=2.25,  t  =0.09  Inch,  a  =14  KSI, 

M.1,  5d-008  Inch  °  1  B 
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